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(57) ABSTRACT 

A method for corrections of measurements of points of inter- 
ests measured by beams of radiation propagating through 
stratified media including performance of ray-tracing of at 
least one ray lunched from a metrology instrument in a direc- 
tion of an apparent point of interest, calculation a path length 
of the ray through stratified medium, and determination of 
coordinates of true position of the point interest using the at 
least one path length and the direction of propagation of the 
ray. 
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LASER RADAR THROUGH THE WINDOW 
(LRTW) COORDINATE CORRECTION 
METHOD 

STATEMENT REGARDING FEDERALLY 5 

SPONSORED RESEARCH AND DEVELOPMENT 

Parts of this invention were conceptualized developed and 
reduced to practice using US government funding under 
National NASA contract for MSES II/A project #07530 10 
“JAMES WEBB SPACE TELESCOPE (JWST) INTE- 
GRATED SCIENCE INSTRUMENTAL MODULE (IS- 
TIM) METROLOGY”. The US government has certain 
rights in this invention. 15 

RELATED APPLICATIONS 

The current application results from and claims priority 
benefit of U.S. Provisional Application Ser. No. 61/674,985, 2 o 
entitled “LASER RADAR THROUGH THE WINDOW 
(LRTW) COORDINATE CORRECTION METHOD”, filed 
Jul. 24, 2012, which is herein incorporated by reference in its 
entirety. 

25 

FIELD OF THE INVENTION 

The invention relates generally to methods for correction 
of physical effects imparted on a beam after it has been 
transmitted through different media (also referred to as envi- 30 
romnents or interfaces) to measure targets of interest. More 
particularly, the present invention pertains to measurements 
by propagation of radiation beams through stratified media 
(media that may be modeled via superposition of multiple 
layers ofmaterials having substantially similar physical prop- 35 
erties pertinent to individual layers while the physical prop- 
erties my vary across the layer boundaries (e.g. when a beam 
of radiation crosses over the boundaries between the stratified 
layers). Even more particularly, several embodiments of the 
present invention pertain to optical measurements through 40 
transparent windows of vacuum chambers or vehicles (either 
manned or unmanned) such as consisting of satellites, rock- 
ets, missiles, atmospheric vehicles, exoatmospheric vehicles, 
aircrafts, spacecrafts, surface vehicles, underwater vehicles, 
combination vehicles (e.g. “amphibians”), and similar. 45 

BACKGROUND OF THE INVENTION 

It should be noted that, at least regarding the current inven- 
tion, metrology instruments include measuring instruments 50 
arranged to determine coordinates or relative positions of 
predetermined points of interest using photonic beams and 
rays (visible, infrared, ultraviolet, x rays, gamma rays, micro- 
waves, radio waves), particulate beams (including electron 
and positron beams, neutron beams, pion (charged and neu- 55 
tral) beams, meson (charged and neutral) beams, ion beams, 
atomic and polyatomic beams, molecule beams, and combi- 
nations including neutral and charged plasma beams and jets, 
and/or beams and streams of aggregated structures (including 
but not limited to, colloidal complexes, mono and poly crys- 60 
tals, nanotubes, fullerenes, glassy aggregates and combina- 
tion structures). 

Consequently, a particular subset of optical metrology 
instruments may include Laser Radars (LR), Laser Trackers 
(LT), optical directional and rangefinders, optical imagers 65 
and cameras, optical scanners, laser levels and plane finders, 
goniometers, theodolites, clinometers, tiltmeters, optical 


2 

sights, optical markers and designators, and/or combination 
metrology instruments combining structures and functional- 
ities of the listed instruments. 

In particular embodiments, optical metrology instruments 
such as LR or LT may have measurement accuracies on the 
order of 0.025 nun over a 2 m measurement range when used 
in ambient conditions with no window between the measure- 
ment instrument and the target. When these same instruments 
are used to measure taigets through a thick dielectric window 
without correcting for window effects measurement error will 
be on the order of the thickness of the single-pass glass path. 
Consequently, measurements made through a window, with- 
out applying the necessary corrections may result in subop- 
timal measurements. 

In different classes of embodiments addressing situations 
where the metrology instrument and object to be measured 
may be separated by a stratified medium because of environ- 
mental, safety, contamination, or packaging requirements, 
etc. could potentially use the LRTW correction capability. 
One such example may consider accurate remote measure- 
ments and recordings of positions of parts (e.g. fuel roads) 
inside nuclear reactors trough windows, cooling fluids, and/ 
or moderating fluids. 

Prior to the development of the LRTW coordinate correc- 
tion method, correction of LRTW measurements have been 
performed generally using connnercial ray trace methods and 
software usually augmented with custom scripts. Frequently, 
however, making through-the-window measurements with 
metrology instruments such as the LR or LT was avoided 
altogether. 

Commercial ray trace methods and software packages may 
be costly and post -processing of the ray trace results may still 
be required. Also, calling commercial ray trace applications 
from a commercial metrology application, although possible, 
may involve additional overhead for importing/exporting the 
data and for running uncompiled scripts. 

The LRTW correction software of the current invention 
may be integrated into connnercial metrology software to 
allow targets measured through a window to be corrected in 
seconds to give the user real-time feedback. This LRTW 
correction software offers a complete solution for correcting 
all range and pointing errors associated with measuring tar- 
gets through a window. Also, costly connnercial ray trace 
software may not be required since this functionality may be 
an integral function of the LTRW software. 

The LRTW coordinate correction method may also be used 
in a reverse mode to calculate apparent point of interest coor- 
dinates for known or corrected taigets measured through a 
window. The apparent point of interest coordinates refers to 
the instrument reported coordinates when no window, 
vacuum, or SGR corrections have been applied. This capabil- 
ity may be useful for allowing the metrology instrument to 
locate targets measured through a window based on ambient 
blue print values to within 0.050 mm. Without this capability 
the time required to locate the targets manually could increase 
substantially as targets inside a vacuum chamber may not be 
well illuminated and may be difficult to locate, especially 
when viewing through the built-in LR or LT camera. 

Some aspects and issues related to the subject matter of 
current invention have been addressed in the following non- 
patent publications, incorporated herein by reference: 

T. Hadjimichael et. al. in “Cryogenic metrology for the 
James Webb Space Telescope Integrated Science Instrument 
Module alignment target fixtures using laser radar through a 
chamber window”, Proc of SPIE Vol. 7793. 

B. Eegholm, T. Hadjimichael, J. Hayden. R. Ohl, D. 
Kubalak, R. Telfer, “Laser Radar Through the Window 
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(LRTW) Coordinate Correction Software”, Sigma Space 
Invention Disclosure, filed Nov. 4, 2010. 

SUMMARY OF THE INVENTION 

Laser Radar measurements may be made through a 
vacuum chamber window for cryogenic metrology of targets 
attached to two different alignment fixtures used for the 
James Webb Space Telescope. The range and pointing errors 
imparted on the LR beam transmitted through up to four 
different environments resulted in reported LR distance 
errors of up to 22 mm, or roughly half the total single-pass 
glass path. These environments included: the air outside the 
chamber where the LR resides, the chamber window, the 
chamber environment, and the glass forming a solid glass 
retroreflector (SGR) target. The custom correction method 
and corresponding MATLAB® software have been devel- 
oped to correct the measured laser-radar-through-the-win- 
dow coordinates in post-processing. The software employs 
real ray tracing combined with a closed-form optical path 
length solution to resolve pointing and ranging errors caused 
by the phase- and group -refractive indices, respectively. The 
optical prescription used for ray tracing is established with 
knowledge of the environmental and physical parameters 
associated with each of the environments. The LRTW correc- 
tion code shows agreement with independently developed 
FRED optical ray trace scripts to better than 1 pm average 
distance for identical environmental input parameters. The 
LRTW correction code was also verified by comparing 
LRTW corrected targets at ambient temperature and pressure 
on the JWST alignment fixtures to those same targets when 
measured with no window in place. The resulting average 
distance error was less than 30 pm (2-sigma), which is con- 
sistent with LR-related sources of measurement error for this 
test configuration. The LRTW coordinate correction software 
could also be used “in reverse” to calculate apparent point of 
interest locations for LR measurements made through a win- 
dow given the blueprint or nominal values of those targets. 
This capability is useful for allowing the LR to quickly locate 
targets measured through a window or when setting up scan 
boundaries around poorly illuminated distant targets. 
Although the Laser Radar (LR) is referenced throughout this 
document, it is implicit that measurements made through a 
window by any time-of-flight or line-of-siglit based optical 
metrology instrument could be corrected using the LRTW 
code. 

The present invention is also directed to devices and assem- 
blies arranged to perform a set of coordinate correction steps 
(or any individual step as required) and to perform error 
analysis and error assessment and management appropriate 
for different implementations of the method of current inven- 
tion in various embodiments. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. lisa schematic of a test configuration for measuring 
the LR point of interest inside a vacuum chamber with the LR 
in accordance with one embodiment of current invention. 

FIG. 2 is a schematic of the direct & through measurements 
used on the first surface of the uncoated chamber window to 
establish the position and normal of the window for the ray 
trace, where the LR, at P0, measures a real tooling ball at PI 
and measures the virtual TB at point P2. 

FIG. 3 is a plot of microstrain versus temperature for Schott 
N-BK7 glass determined from average of ESA data (Excel 
file name “BK7 at PMIC low temp.xls”, ESA contract 1358, 
from testing on 6 Nov. 2009 and 10 Nov. 2009). 


4 

FIG. 4 is an example of plots of absolute phase and group 
refractive indices versus temperature for a) N-BK7, and for b) 
Coming Fused Silica 7980, generated using the wavelength- 
and temperature-dependent Sellmeier equation with corre- 
5 sponding CHARMS coefficients. 

FIG. 5 is a schematic of a test configuration for LR mea- 
surements in accordance with another embodiment of current 
invention. 

FIG. 6 is a schematic of a test configuration for LR mea- 
1 0 surements in accordance with another embodiment of current 
invention. 

FIG. 7 is a schematic of a test configuration for LR mea- 
surements in accordance with another embodiment of current 
15 invention. 

DETAILED DESCRIPTION 

In the following description of embodiments of the present 
, 0 invention, numerous specific exemplary details are set forth 
in order to provide a more thorough understanding of the 
present invention. However, it will be apparent to one skilled 
in the art that the present invention may be practiced without 
one or more of these exemplary details. In other instances, 
25 well-known features of prior art have not been described in 
detail to avoid unnecessarily complicating the description. 

One embodiment of the current invention utilizes a MAT- 
LABs code developed to model and correct for the refraction 
and range effects associated with making LR measurements 
30 through multiple optical media. Real ray tracing was used to 
model the refractive effects on the LR beam using known or 
calculated information about each of the media. In particular 
embodiments, a manually populated Excel spreadsheet was 
used to supply the LRTW correction software with the 
required optical, environmental, material, and dimensional 
parameters required to correct the apparent xyz-points mea- 
sured by the LR. The LRTW correction software generates a 
unique optical prescription for each apparent point of interest 
40 to be corrected. The optical prescription, containing i=5 sur- 
faces, is in the general form Z t , R„ n,, and, N i; where Z is the 
global position along the Z-axis at which the surface inter- 
sects the Z-axis, R is the radius of curvature, n is the phase 
refractive index, and N ; is the unit normal directions of the 
45 SGR entrance face. 

In this embodiment, a fixed-position LR mounted outside a 
vacuum chamber window was positioned with its nodal point 
at the origin of a coordinate system, called the Ray Trace 
Frame, whereby the Z-axis was defined to be normal to the 
50 air-side of the chamber window and the Y-axis is pointing in 
the gravity direction. The LR was roughly centered on the 
window so that it had a line-of-sight through the window to 
the targets inside the chamber. An LR beam with a wave- 
length of 1550 mn was transmitted through several different 
55 transparent optical media while measuring multiple targets 
inside the chamber. The transparent optical media included 
(but have not been limited to): 1 ) the ambient air around the 
LR outside the chamber, 2) the chamber window, and 3) the 
ambient or cryogenic vacuum chamber environment. It may 
60 be useful to consider additional (e.g. fourth) optical media, 
for example in cases where the point of interest may include 
the apex of a solid glass retroreflector (SGR) made of N-BK7. 
Real ray tracing may be used to trace a ray from the LR nodal 
point through each of the optical media to the true point of 
65 interest position V T . 

.An exemplary ray trace configuration 1 00 is shown in FIG. 
1 . For completeness, the SGR 105 is shown in this ray trace 
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configuration 100; however, if the point of interest is not an 
SGR 105 then the SGR 105 thickness may be set to h=0 so 
that V r =P 3 . 

Hie custom MATLAB® code [MATLAB® 64-bit version 
8.1.0.604 (R2013a) by MathWorks®] of this embodiment 5 
uses well established vector-based refractive ray tracing 
methods [e.g. Foley, J. D., and Van Dam, A. “Computer 
Graphics: Principles and Practice, Second Edition” Addison- 
Wesley New York, N.Y. 1990.]; Glassner, A. (ed) “An Intro- 
duction to Ray Tracing” Academic Press New York, N.Y. to 
(1989)] to launch rays from the LR through a simple optical 
system. The ray trace may be initiated by calculating the 
direction of the ray exiting the nodal point 110 as L, - 

ULi/UiA! where the LR nodal point 110 is P o =<0, 0, 0> 
and the apparent point of interest 115 in the Ray Trace Frame 15 
is V A =<V Ax , V Ay , V Az >- The apparent point of interest 115 
position may be the position of the point of interest as mea- 
sured and reported by the LR when the refraction and range 
effects associated with beam propagation through multiple 
optical media are not considered. 

The true point of interest 120 position, V T , cannot be 
resolved with ray tracing alone. However, by temporarily 

setting Z 3 =V 4z , it is possible to resolve L, L, L 3 , and Z 4 . 

Note that only the direction of L } is known so V r cannot be 
determined until La is calculated with knowledge of the cor- 
rected range. 

The distance Z l between the LR, located at P0, to the first 
surface of the window 125,Sl,canbe determined using direct 
and through measurements [Spatial .Analyzer Users Manual, 30 
New River Kinematics, Williamsburg, v. 1.21.2008, page 
150,] FIG. 2, on a single tooling ball (TB) 205, placed at PI, 
and its refection 210 defined by the virtual point P2. The TB 
205 position, although somewhat arbitrary, must satisfy the 
minimum LR to TB distance on the order of 1 m and the LR 
must have lines of sight to both PI and P2. The measurement 
of the virtual TB reflection 210 position at P2 requires that the 
signal detection parameters in the LR control software be 
modified to accept a low reflectance target. The low reflec- 
tivity from P2 results from the LR beam double-passing the 411 
uncoated 4% reflective window. The window surface 125, S 1 , 
is normal to the Z axis by definition of the ray trace frame. 

The unit vector, n, describing the normal of the first surface 
125, SI, of the window is 


PI -PI (1) 

" = |P1 - Pl \ ' 

50 

where PI and P2 are the real and virtual tooling ball (TB) 
points, respectively, measured by the LR. 

The vacuum chamber window 130 will undergo an axial 
deflection as a function of differential pressure that imparts a 
meniscus shape on the window 130, FIG. 1, where the radius 55 
of curvature, R, for either side of the window is described by 
Hearn [D. Hearn, “Vacuum window optical power induced by 
temperature gradients”, SPIE Vol. 3750, July (1999)] as 

R=(4Eh 3 )/[3(l-v)(3+v)a 2 Ap] (2) gQ 

where G, is the elastic modulus, h is the window thickness, v 
is Poisson’s Ratio, a is the radius of the simply supported 
window, and Ap is the pressure differential on the window 
130. In one embodiment, a 25.502 nun thick Fused Silica 
window 130 used in a vacuum chamber 135 may result in 65 
approximately 400 m radius of curvature for a differential 
window pressure of Ap=10 -5 mm of Hg. The radius of cur- 


vature of the window 130 may be calculated, using equation 
(2), for each apparent point measured by the LR based on the 
differential pressure on the window 130 at the time of the 
measurement. 

In the exemplary embodiment of FIG. 1, using a tilted 
plane-parallel plate model of the SGR 105, may be used by 
requiring that all rays entering the SGR 105 end at the apex, 
or mechanical vertex, of the SGR 105. Some motivation for 
such constraint may be found by considering that the SGR 
105 apex may be intrinsically identified when measured with 
a laser tracker LT, which automatically locks-on to the SGR 
105 apex. Consequently, only the chief ray of the LT beam 
entering the SGR 105 needs to be raytraced. In the case of the 
LR, identification of the SGR 105 apex may be complicated 
by the need to use post-processing to establish the apex from 
the SGR 105 facet intersections. Surface normal 140 of the 
SGR entrance face may be determined from separate inde- 
pendent measurements on the adjacent co-planer mount sur- 
face of the SGR 105 that are assumed to be the same at 
ambient and cryogenic temperature. The mechanical vertex 
may be used as the reference for the SGR 1 05, rather than the 
SGR 105 nodal point, to facilitate comparisons with direct 
back-side measurements of the SGR 105 mechanical vertex 
from the 3-plane intersection of the facets. 

A thickness versus temperature correction can be also 
applied to the SGRs using empirical data from the European 
Space Agency (ESA) where microstrain was defined as a 
function of temperature, FIG. 3, using the polynomial 



(3) 


where P 0 =- 148.53, Pj = 6 . 7335333, P 2 =0.0014505, andP 3 =- 
2.821 5e-005. Using this calculation our 1 9 mm thick N-BK7 
SGRs change thickness by -24 pm as the temperature 
changes from 20 C to -238 C for an initial thickness of 1 9 mm 
at 20° C. 

High-accuracy time-of-flight instrumentation such as the 
LR and LT may use Edlen’s equation [K. R Birch and M. J. 
Downs, “An updated Edlen equation for the refractive index 
of air", Metrologia 30, 155-162 (1993) K. P. Birch and M. J. 
Downs, “An updated Edlen equation for the refractive index 
of air,” Metrologia 30, 1 55-1 62 (1 993); K. P. Birch and M. J. 
Downs, “Correction to the updated Edlen equation for the 
refractive index of air”, Metrologia 31, 315-316 (1994); G. 
Bonsch and E. Potulski, “Measurement of the refractive index 
of air and comparison with modified Edlen’s formulae”, 
Metrologia 35, 133-139 (1998)], to calculate the absolute 
phase refractive index of air as a function of the air pressure, 
temperature, and relative humidity. The LR may be equipped 
with an external environmental sensor that provides pressure, 
temperature, and relative humidity input to the Edlen’s equa- 
tion for all measurements made by the LR. However, in the 
case where the vacuum chamber 135 contains air which is 
separated from the air 145 outside the chamber, a separate 
phase refractive index may be calculated using Edlen’s equa- 
tion using input from the appropriate sensors placed inside 
the air-filled vacuum chamber. The phase refractive index 
may be used with the real ray tracing to model refractive 
effects at the interfaces. 

The absolute phase refractive index of the glass window 
130 and the glass SGR 105 maybe described by the tempera- 
ture- and wavelength-dependent modified Sellmeier disper- 
sion equation [D. Leviton, B. Frey, “Temperature-dependent 



US 9,109,898 B2 


7 


8 


absolute refractive index measurements of synthetic fused 
silica”, Proc of SPIE Vol. 6273; D. Leviton, B. 1. Frey, “Cryo- 
genic, High-Accuracy, Refraction Measuring System 
(CHARMS) — a new facility for cryogenic infrared through 
far-ultraviolet refractive index measurements”, SPIE, Glas- 
gow, June, 2004], using 


n( A, T) = 



S,(T )-\ 2 
A 2 -Xf(T) 


(4) 


where 


-continued 

where 
k = 1 /A. 

5 

The group index for the air may be calculated using a 
closed-form solution for dn/dk, from Edlen’s equation, to 
model range effects. The group refractive indices for the 
I0 Coming Fused Silica 7980 chamber window and the N-BK7 
SGR are calculated from the derivative of the phase refractive 
index in equation (4) to get 


4 

S i (T) = Y j S u -T J , 

7=0 

and 


(5) 

15 


3 

dn k SjXj 

dk n / i (\-A.:k 2 ) 2 ' 

j - 1 


( 8 ) 


4 (6) 

A;(7') = 2> ir 7T 
;=o 

Determination of the true point of interest 120 position, 
in FIG. 1, from the LR-measured apparent point of interest 
115 position, V A , requires that real ray tracing be combined 
with a closed-form range correction equation. The range cor- 


The absolute phase refractive index, n, and absolute group 
refractive index, n g , FIG. 4, of the Corning Fused Silica 7980 
chamber window and the N-BK7 SGR are calculated using 
equations (4)-(8) at the fixed LR wavelength of 1550 mil. 
equations (4), (5), (6), and (8) are based on the modified 
Sellmeier equation and the CHARMS coefficients in Table 1 
(a). 


TABLE 1 


Coefficients for the temperature-dependent 
Sellmeier fit of the refractive index of: a) Coming Fused Silica 7980, 
valid for 30 K < T < 300 K, and 0.4 pm < A, < 2.6 pm, and b) N-BK7, valid 
for 50 K < T < 300 K. and 0.45 inn < A. < 2.1 pm. 


Term 

SI 

S2 

S3 

LI 

L2 

L3 

(a) CHARMS coefficients for Coming Fused Silica 7980 glass 

-j-0 

1.10127E+00 

1.78752E-05 

7.93552E-01 

8.90600E-02 

2.97562E-01 

9.344S4E+00 

T 1 

-4.9425 IE-05 

4.76391E-05 

-1.27815E-03 

9.08730E-06 

-8.59578E-04 

-7.09788E-03 

T 2 

5.27414E-07 

-4.490 19E-07 

1.845 95 E-05 

-6.53638E-08 

6.59069E-06 

1.01968E-04 

T 3 

-1.59700E-09 

1.44546E-09 

-9.20275E-08 

7.77072E-1 1 

-1.09482E-08 

-5.07660E-07 

-j-i 

1.75949E-12 

-1.57223E-12 

1.48829E-10 

6. 84605 E- 14 

7.85145E-13 

8.21348E-10 



(b) CHARMS coefficients for N-BK7 (or BK7) glass 


•pO 

1.0703 6E+00 

2.02896E-01 

8.50659E-01 

8.20884E-02 

1.41222E-01 

9.38638E+00 

T 1 

1.75658E-03 

-1.75289E-03 

1.03233E-03 

4.35600E-05 

-5.41379E-05 

5.245 15E-03 

T 2 

-4.28401E-06 

4.18934E-06 

-6.88935E-06 

5.85421E-08 

4.72890E-06 

-3.54202E-05 

T 3 

-1.65687E-11 

4.52824E-10 

2.565 10E-08 

-4.13612E-10 

-2.73546E-08 

1.32249E-07 

-p4 

7.55467E-12 

-8.18792E-12 

-4.71564E-11 

5.61059E-13 

8.01380E-11 

-2.38407E-10 


rection equation uses the results of the real ray trace, namely 

i, , L, L 3 , and L t , combined with the calculated group 50 
refractive indices for each optical media, to solve for the true 
point of interest position. 

When a pulse of light from either an LR or LT is propagated 
through a dispersive media the modulation envelope 
advances at a velocity that is different from the phase velocity 55 
of the waveform. Consequently, the group refractive index 
[“Optics — A rh Edition”, Eugene Hecht, p 296, Ch 7.2.2 Group 
Velocity”] may be used when modeling the range effects. 
Both the phase- and group-index of light propagating through 
a dispersive medium have a wavelength and temperature 60 
dependence. For a homogeneous dispersive media the group 
index n g is related to the phase index n, by 


The true total geometric path length of the LR beam as it 
travels from P 0 to V r through four optical media is 

' /, | +/, ^+/, ,+/,4 (9) 

where L ( , for i=l , 2,3,4, is the geometric length of the ray 
in each of the optical media. The optical path length measured 
by the LR is 

L A ng 1 =L 1 Kg 1 +Z, 2 Kg 2 +Z, 3 Kg 3 +L 4 Kg 4 , (10) 

where L A is the apparent total geometric length of the LR 
beam, and ng„ for i=l, 2, 3, 4 is the group index of each of the 
optical media. Solving for L 3 gives the closed-form range 
equation 


(L A ng l - L l ng l - hing 2 - Ung 4 ) 


(7) 65 


(ID 
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Finally, using equation (11) with the values L,, L 2 , L 3 , and 

L 4 determined from the ray trace, and using the calculated 
group refractive indices for each of the optical media, the true 
point of interest coordinates (V T x, V r y, V jZ) can be resolved 
by vector addition of the rays in each of the optical media as 

ty= L , + L 2 +Z 3 Z. 3 + 4, . ( 12 ) 

Some additional functionalities of different embodiments 
of the current invention include calculation of apparent point 
of interest 115, V^, given the true or nominal point of interest 
positions 120, V r (having known coordinates (V r X, V r y, 
VjZ)). This capability, FIG. 5, may be useful for allowing the 
LR to efficiently locate measured through a stratified optical 
medium (e.g. window 130) or when setting up scan bound- 
aries around poorly illuminated distant targets. Occasionally, 
it may be also helpful to convert LR measurements made on 
metrology targets in a lab to their corresponding values in a 
particular vacuum chamber environment when measured 
through a window. 

The apparent point of interest 115 position, V A , can be 
calculated by using the LRT W code in an optimization loop to 
minimize the merit function 

/= Wj-Vr’I. ( 13 ) 

where V r ' is the current calculated point of interest 520 posi- 
tion using guess V A ' for the apparent point of interest 515 
position, orV r '=LRTW(V^'). The optimization loop modifies 
the apparent point of interest 515 position V A until the merit 
function has been optimized (e.g. until f„ ! „>IV 7 ^V 7 .'l, 
wherein the f mi „ represents a predetermined optimization 
parameter corresponding to the desired accuracy). The opti- 
mization may be initialized with V^-Vj. as V r may be suffi- 
ciently close to V A to allow rapid convergence. 

The above methods for determination of the true point of 
interest 120 position, (as illustrated, for example, in FIG. 1.) 
and determination of the apparent point of interest 115 posi- 
tion, V A , (as illustrated for example, in FIG. 5.) may be 
schematically represented in flowcharts 600 find 700 forms in 
respective FIGS. 6 and 7. Regarding the flowchart 600 rep- 
resenting the exemplary embodiment of the method for deter- 
mination of the true point of interest 120 position, initial steps 
of determination of components and physical parameters 610 
and definition of the set-up geometry 620 may be initially 
performed during the conceptualization and set-up phases. 
One should note that steps 610 and 620, depending on par- 
ticulars of different embodiments may be performed in 
sequence (in direct or reversed order) or iteratively (e.g. on 
the component by component basis) during the conceptual- 
ization and definition of the model constituencies and geom- 
etry. It should be also evident to practitioners (at least from 
disclosures above) that increasingly detailed description of 
the optical and special relationships of the relevant param- 
eters (e.g. detailed geometry and properties of the chamber 
window 130 and surface 125 under various pressures and 
temperatures) may be introduced and/or subsequently refined 
without exceeding the scope and limitations of the current 
invention. 

Regarding the performance of the ray-tracing step 630, one 
may note that, depending upon particular embodiment, dedi- 
cated ray -tracing routines may be utilized, combining com- 
mercially available and/or custom developed programs and/ 
or program modules. 

Regarding the determination of coordinates steps 650 of 
various embodiments, one may note that the disclosed 
method may be repeated in its entirety or in parts (e.g. steps 
630-650) such that a plurality of points of interests may be 


obtained and used for generation of images and 3D object 
positions, for example for imaging, images corrections, 3D 
orientation, guidance, position and motion control, mapping, 
reconnaissance, and other embodiment-dependent usages. 

5 Regarding the embodiments of the method for determina- 
tion of the (unknown) apparent point of interest 115 position, 
V A , corresponding to the (given, known or previously deter- 
mined for example as in step 710) true point of interest 120 
position, V,) as schematically represented in the flowchart in 
to FIG. 7, it may be noted that the embodiment-specific steps 
610 and 620 may not differ substantially for the analogues 
steps as represented in FIG. 6. 

Regarding the step 720 of determination of the current 
iteration calculated point of interest 520 position V r ' it may be 
15 noted that the aforementioned initial choice of iteration start- 
ing point V_/=V r of the current calculated point of interest 
520 position V r ' may be convenient in some embodiments, at 
least in part, because of the reasonable expectation of prox- 
imity of Vj and V r . In different embodiments, the iterative 
20 process represented by the flowchart 700 may be initiated by 
a semi-arbitrary or other convenient guess for the apparent 
point of interest 515 position V 4 '. Correspondingly, the step 
750 of definition on new updated choice of apparent point of 
interest V A ’ may be conducted using various embodiment- 
25 specific methodologies generally motivated by the practitio- 
ners’ interest to accelerate and improve the iteration process 
(e.g. to efficiently and reliably satisfy the convergence crite- 
rion exemplified by the interrogation block 740). 

Also, the step 760 of determination of the apparent point of 
30 interest 115 position from the convergence of the iterative 
process may be followed by subsequent repetitions of the step 
710 for an additional (e.g. proximal) point of interest 515 such 
that subsequent iterative procedures may be executed for a 
plurality of additional points of interest 515. 

35 Targets such as pinholes, thru-holes, and SGRs used on 
JWST alignment fixtures could be successfully measured 
through a window and corrected for range and pointing 
errors. The LRTW correction code was verified by comparing 
the LRTW corrected targets on the JWST alignment fixture 
40 targets with the same targets measured with no window in 
place. The resulting average distance error was less than 30 
pm (2-sigma), which is consistent with LR-related sources of 
measurement error for this test configuration. The LRTW 
correction code showed agreement with independently devel- 
45 oped FRED optical ray trace scripts to less than 1 pm average 
distance for identical environmental input parameters. The 
LRTW coordinate correction software can also be used “in 
reverse” to calculate apparent point of interest locations for 
LR measurements made through a window given the blue- 
50 print or nominal values of those targets. The LRTW correc- 
tion capability has been reduced to practice using both laser 
radar and laser tracker measurement devices. Although the 
LRTW correction was initially developed to correct for laser 
radar and laser tracker measurements made through a win- 
55 dow. this correction technology could be applied to any 
metrology device that uses pointing and/or ranging to mea- 
sure the position of a metrology point of interest through a 
window. 

The present invention has been described with references 
60 to the above exemplary embodiments. While specific values, 
relationships, materials and steps have been set forth for 
purpose of describing concepts of the invention, it may be 
appreciated by persons skilled in the art that numerous varia- 
tions and/or modifications may be made to the invention as 
65 shown in the specific embodiments without departing from 
the spirit or scope of the basic concepts and operating prin- 
ciples of the invention as broadly described. It should be 
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recognized that, in the light of the above teachings, those 
skilled in the art can modify those specifics without departing 
from the invention taught herein. Having now fully set forth 
the preferred embodiments and certain modifications of the 
concept underlying the present invention, various other 
embodiments as well as certain variations and modifications 
of the embodiments herein shown and described will obvi- 
ously occur to those skilled in the art upon becoming familiar 
with such underlying concept. It is intended to include all 
such modifications, alternatives and other embodiments inso- 
far as they come within the scope of the appended claims or 
equivalents thereof. It should be understood, therefore, that 
the invention may be practiced otherwise than as specifically 
set forth herein. Consequently, the present embodiments are 
to be considered in all respects as illustrative and not restric- 
tive. 

The scope of the LRTW capability is not limited to optical 
metrology measurements on telescope optics or hardware 
made through a chamber window. Further improvements 
could be made to the implementation of LRTW corrections 
by employing multiple metrology devices used through a 
window in a unified spatial metrology network (USMN) to 
reduce measurement errors. 

We claim: 

1. A method for corrections of metrology instruments mea- 
surements of positions of points of interest measured by 
beams of radiation propagating through stratified media com- 
prising following steps: 

a) definition and arrangement of at least one metrology 
instrument, definition of at least one point of interest, 
and at least one stratified medium; 

b) definition and arrangement of at least one measurement 
setup having the at least one metrology instrument, at 
least one point of interest, ad at least one stratified 
medium having a least a portion of stratified medium 
arranged in between and separating the at least one 
metrology instrument and the at least one point of inter- 
est; 

c) performance of ray-tracing of at least one ray lunched 
from that at least one metrology instalment in a direction 
of an apparent position of the at least one point of interest 
through at least the portion of the at least one stratified 
medium and determination of at least one direction of 
propagation of the at least one ray through the at least the 
portion of the at least one stratified medium arranged in 
between and separating the at least one metrology 
instrument and the at least one point of interest, and 
measurement of at least one apparent total geometric 
length from the at least one metrology instalment to the 
apparent position of the at least one point of interest 
along the at least one ray; 

d) calculation of at least one path length of the at least one 
ray through the at least a portion of the at least one 
stratified medium using optical properties of the at least 
one stratified medium, the at least one direction of 
propagation of the at least one ray, and the at least one 
apparent total geometric length from the at least one 
metrology instalment to the apparent position of the at 
least one point of interest along the at least one ray; and 

e) determination of coordinates of taie position of the at 
least one point interest using the at least one path length 
and the at least one direction of propagation of the at 
least one ray. 

2. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim. 1, 
wherein the at least one metrology instalment has been 
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arranged to radiate a beam of radiation through the at least a 
poaion of the: at least one stratified medium. 

3. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 2, 

5 wherein the at least one metrology instalment has been 
arranged to radiate a beam of radiation chosen from a set of 
beams consisting of visible photon beams, infrared, photon 
beams ultraviolet photon beams, x-ray photon beams, gamma 
ray photon beams, microwave, photon beams, radio waves 

10 photon beams, particulate beams, electron, beams, positron 
beams, pion beams, meson beams, proton neutron beams, 
alpha particle beams, ion beams, atom beams, polyatomic 
beams, molecule beams, neutral and changed plasma beams, 

15 colloidal complexes beams, mono and poly crystals beams, 
nanotubes beams, fullerenes beams, glassy aggregates beams 
and combinations of listed beams. 

4. The method for corrections of metrology instruments 
measurements or positions of points of interests of claim 1 , 

20 wherein the at least one metrology instrument is an optical 
metrology instrument chosen from the set of optical metrol- 
ogy instruments consisting of laser radars, laser trackers, 
optical rangefinders, optical imagers and cameras, optical 
scanners, laser levels and plane finders, goniometers, theodo- 

25 lites, clinometers, tiltmeters, optical sights, optical markers 
and designators, and combination metrology instruments 
combining structures and functionalities of the listed optical 
metrology instruments. 

5. The method for corrections of metrology instruments 

30 measurement of positions of points of interests of claim 1 , 

wherein the at least one stratified medium is a dielectric 
window. 

6 . The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 5, 

35 wherein the dielectric window have been arranged on a 
vacuum chamber. 

7. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 5, 
wherein, the dielectric window have been arranged on a 

40 vehicle. 

8. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 7, 
wherein the vehicle is chosen from a set of vehicles consisting 
of satellites, rockets, missiles, atmospheric vehicles, exoat- 

45 mospherio vehicles, aircrafts, spacecrafts, surface vehicles, 
underwater vehicles, and combinations. 

9. The method for corrections of metrology instruments 
measurements of positions of point of interests of claim 1, 
wherein the at least one stratified medium includes at least 

50 one solid glass retroreflector (SGR) target. 

10. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 9, 
wherein the at least one glass retroreflector (SGR) target has 
been arranged in proximity to the at least one point of interest. 

55 11. The method for corrections of metrology instruments 

measurements of positions of points of interests of claim 1 
wherein after the step e) of determination of coordinates of 
taie position of the at least one point of interest, the steps c)-e) 
are repeated in succession for determination of coordinates of 

60 taie position of the at least one point of interest. 

12. A method for corrections of metrology instalments 
measurements of positions of points of interests measured by 
beams of radiation propagating through stratified media com- 
prising following steps: 

65 a) definition and arrangement of at least one metrology 
instrument, definition of at least one point of interest, 
and at least one stratified medium; 
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b) definition and arrangement of at least one measurement 
setup having the at least one metrology instrument, at 
least one point of interest, and at least one stratified 
medium having at least a portion of stratified medium 
arranged in between and separating the at least one 5 
metrology instrument and the at least one point of inter- 
est; 

c) determination of coordinate of a true point of interest 
position of the at least one point of interest; 

d) performance of ray-tracing of at least one ray lunched 10 
from that at least, one metrology instrument in a direc- 
tion of a choice point of interest position through at least 
the portion of the at least one stratified medium and 
determination of at least one direction of propagation of 
the at least one ray through the at least the portion of the 15 
at least one stratified medium arranged in between and 
separating the at least one metrology instrument and the 

at least one point of interest, and measurement of at least 
one apparent total geometric length from the at least one 
metrology instrument to the apparent position of the at 20 
least one point of interest along the last one ray; 

e) calculation of at least one path length of the at least one 
ray through the at least a portion of the at least one 
stratified medium using optical properties of the at least 
one stratified medium, the at least one direction of 25 
propagation of the at least one ray, and the at least one 
apparent total geometric length from the at least one 
metrology instrument to the apparent position of the at 
least one point of interest along the at least one ray; 

f) determination of coordinates of a current calculated 30 
point of interest position using the at least one path 
length and the at least one direction of propagation of the 

at least one ray; 

g) determination of at least one merit function value by 
comparing separation of the current calculated point of 35 
interest position and the true point of interest position. 

h) coordinate of the at least one merit function value at least 
one predetermined optimization parameter if the at lean 
one predetermined optimization parameter substantially 
exceeds the at least one merit function value the coordi- 40 
nates of the current calculated point of interest position 
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are found to represent the coordinates of an apparent 
point of interest position corresponding to the true point 
of interest position the at least one point of interest, and 
if the at least one predetermined optimization parameter 
does not substantially exceed the at least one merit func- 
tion value a choice point of interest position is updated 
and the step d)-h) repeated in sequence. 

13. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 12, 
wherein the at least one stratified medium is a dielectric 
window. 

14. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 13, 
wherein the dielectric window has been arranged on a 
vacuum chamber. 

15. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 13, 
wherein the dielectric window lias been arranged on a vehicle. 

16. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 15, 
wherein the vehicle is chosen from a set of vehicles consisting 
of consisting of satellites, rockets, missiles, atmospheric 
vehicles, exoatmspheric vehicles, aircrafts, spacecrafts, sur- 
face vehicles, underwater vehicles, and combinations. 

17. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 12, 
wherein the at least one stratified medium includes at least 
one solid glass retroreflectox (SGR) target. 

18. The method fox corrections of metrology instruments 
measurements of positions of points of interests of claim 17, 
wherein the at least one glass retroreflector (SGR) target have 
been arranged in proximity to the at least one point interest. 

19. The method for corrections of metrology instruments 
measurements of positions of points of interests of claim 12, 
wherein after the step e) of determination of coordinates of 
true position of the at least onepoint ofinterest, the steps c)-h) 
are repeated in succession for determination of coordinates 
the apparent point of interest position corresponding to the 
true point of interest position of at least another point of 
interest. 



